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Abstract 

A statistical model of resistor network is proposed to describe a poly- 
mer structure and to simulate the real and imaginary components of its 
ac resistivity. It takes into account the polydispersiveness of the mate- 
rial as well as intrachain and interchain charge transport processes. By 
the application of a transfer matrix technique, it reproduces ac resistivity 
measurements carried out with polyaniline films in different doping de- 
grees and at different temperatures. Our results indicate that interchain 
processes govern the resistivity behavior in the low frequency region while, 
for higher frequencies, intrachain mechanisms are dominant. 
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1 Introduction 

The strong variation of electrical conductivity of conjugated polymer films re- 
sulting from chemical doping, as seen in the insulator-to-metal transitions, was 
first observed in polyacetylene [[!]. Since then, a great amount of experimen- 
tal and theoretical studies has been done to elucidate the mechanisms involved 
in the electronic transport phenomena. Initial efforts were dedicated to un- 
derstand the conduction mechanism along a single chain, with the SSH model 
[Ej, a theoretical one-electron framework based on the Hiickel approach, which 
most successfully described the one-dimensional trans-polyacetylene. Doping 
mechanisms associated with chain defects brought into vogue the concepts of 
solitons and polarons as intermediate levels in the electronic structures of a sin- 
gle conjugated polymeric chain The number of such excitonic defects are 
directly related with the doping level: from weak to moderate doping the con- 
duction mechanism is mainly due to electronic hopping, while extended states 
are reached by strongly doped samples. In macroscopic samples, for which 
extremely large conductivity changes can be observed (circa of 15 orders of 
magnitude), interchain processes must be taken into account owing to interac- 
tions brought about by the molecular packing. Theoretically, such interactions 
may be represented by the interchain exchange integral, which contribute sig- 
nificantly to the final conductivity value. 

The major difficulty in the study of bulk conductivity of polymeric materials 
lies in their highly disordered structure. The carrier transport should contain 
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competing intrachain and interchain processes, and the challenge is not only to 
identify the contribution of each process, but also to quantify them. In this 
context, the alternating conductivity technique is an appropriate tool due to its 
spectroscopic character. 

Among other theoretical treatments, statistical models have also been ap- 
plied to describe polymerization processes B, as well as, to study electrical 
properties of polymers [|[ [| . These models are essentially based on percolation 
concepts and have the advantage of simultaneously treating the phenomena 
involved on both a microscopic and a macroscopic scale. Baughman and Shack- 
lette proposed a random resistor network model which successfully explained the 
conductivity dependence on the conjugated length Q. Andrade et al || inves- 
tigated the influence of percolation disorder on the electrical behavior in the dc 
regime. The distinction between intrachain and interchain mechanisms, both 
equally important to bulk conductivity, nevertheless requires a more detailed 
analysis. 

In this paper we investigate the conductivity properties of weakly doped 
polyaniline in the presence of an alternating electric field, from both the theoret- 
ical and experimental perspectives. Ac conductivity measurements were carried 
out in polyaniline (PANI) films at different temperatures and doping degrees. 
The system is theoretically described by a percolation model where the basic 
ingredients responsible for the conductivity are incorporated. The model also 
takes into account the polydispersiveness of the material, the interchain and 
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the intrachain charge transport. The real and imaginary parts of the resistivity 
were determined using a transfer matrix technique. At low frequencies, inter- 
chain processes are more important and dominate the transport mechanism. On 
the other hand, at high frequencies charge transport should be restricted along 
the polymer chains, so intrachain processes should be dominant. Both regimes 
are well described by our model, which also reproduces the experimental results 
in a remarkable way. 

2 Theoretical Model and Numerical Simulations 

The theoretical model introduced in this paper allow us to focus two important 
aspects of the material: its geometrical structure and its electrical conductivity. 

First, we build a strip by juxtaposing TV square cells, each cell with L x L 
sites (see Fig.l). Here, cells are used to conciliate the polymer size distribution, 
obtained experimentally by the Gel Permeation Chromatography method, with 
the transfer matrix technique developed by Derrida et al. . The use of cells 
was applied here to calculate the electrical conductivity without execeeding the 
capacity of the computer's memory. In its original formulation 0, the tranfer 
matrix acts on a random and ramified structure of resistors. Here, however, only 
a non ramified structure is possible. Therefore, we construct a cell which con- 
tains only linear structures with a Gaussian size distribution and fixed density. 
The information of one cell is stored and transferred, by the matrix method, to 
the next cell. 
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An empty lattice site i of the cell is randomly occupied, representing the 
presence of a monomer. We shall call this site the growing tip. The chain size I is 
taken from a Gaussian distribution probability centered in l with dispersion Al. 
An empty site j, nearest neighbor of the growing tip i, is randomly occupied and 
a resistance R c is attributed to the bond connecting these two monomers (sites). 
The site j becomes now the new growing tip. The process is repeated until the 
chosen polymer size is reached. If the polymer cannot reach the specified length, 
it is discarded. Clearly, a self-avoiding scheme is used to mimic the excluded 
volume interaction. 

After the chain has reached the selected size, a search in all of its monomers 
is performed looking for a pair of nearest neighbor monomers (which can belong 
or not to the same chain), not connected by R c . If such a pair is found then 
these monomers are connected, with probability pd, by a RC circuit. This 
RC circuit is a parallel association where C is the capacitance and Ri is the 
interchain resistance . The idea is to incorporate induced charges and hopping 
(or tunneling) mechanisms. As we shall see, pd is related to the doping level 
and we call it the doping parameter. 

Repetition of the procedure described above built new polymer chains inside 
the cell until some concentration k of the occupied lattice bonds is reached. Let 
N i, be the number of occupied lattice bonds, i. e., N b = number of resistances 
R c + number of impedances Z (the equivalent impedance of the circuit RC) and 
N e b the number of empty lattice bonds (corresponding to infinite resistances) 
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then k = N b/(N b + N e b)- If k — 0.5 the system is at the critical bond per- 
colation threshold of the square lattice. A typical cell is depicted in Fig.l. It 
is considered that a bias voltage is applied on both sides, being L the distance 
between the electrodes. 

To determine the longitudinal resistivity p// of the cell, which is the resis- 
tivity in the direction of the current flow, we follow the tranfer matrix method 
. Compared to the two other main approaches - resolution of Kirchoff's equa- 
tions and node elimination, this technique has many advantages. It consumes 
less CPU time, it avoids troubles with dangling bonds or isolated clusters, and 
the conductivity is calculated exactly. Of course, care should be taken when 
passing from one cell to another: horizontal and vertical bonds on the boundary 
must be stored in order to correctly construct the next cell. We used around 
N ~ 100,000 cells in our simulations. 

Simulated curves p// versus / (the frequency of an ac electric field) are 
presented in Fig. 2. They exhibit two plateaus, i. e., two regions in which the 
resistivity is frequency- independent: one at low frequencies and the other at 
high frequencies. Between these two plateaus the resistivity decays with the 
frequency, obeying approximately the relation p// oc f~ a with a = 1.9. It does 
not agree with the measured value a ~ 1.2, which is the slope of the experimental 
curves of Fig. 3. This discrepancy however, disappears when we consider, as 
we shall see later, the interchain resistance Ri dependence with the frequency 
/. Figures 2a, 2b and 2c correspond to simulations performed at concentration 
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k = 0.5 and they show, respectively, the influence of pd, R c and Ri in the 
P//(f) curves. Ri acts in the low frequency plateau, R c in the high frequency 
and pd causes the whole curve to shift up and down. The interpretation of 
these parameters will be discussed below, when the simulated results and the 
experimental curves carried out in PANI films are compared. Figure 2d was 
obtained at fixed values of pd, R c and Ri but for two different concentrations 
k = 0.5 and k = 0.6. Increasing k simply displaces the whole curve of resistivity 
downward in a frequency independent way. This is expected since, to increase 
the concentration k is equivalent to substitute a number of empty bonds (infinite 
resistances) by occupied bonds (finite resistances). Note that the effects of 
increasing k or pd are very similar - both arc frequency independent and decrease 
the resistivity. As the concentration k does not essentially alter the form of the 
curves, in the rest of this paper we restrict the simulations to be performed at 
k = 0.5. 

3 Experimental Resistivity 

Responses to alternating electrical field applied in doped conjugated polymers 
have been an important subject of research when charge transport phenomena 
arc involved. Some ac results, conductivity versus frequency, have already been 
reported in the literature, with polyacetylene ]8|, £|, poly(3-methylthiophene) 
[ ^0| and poly(o-methoxyaniline) pJj ]. Here we present some results carried out 
with polyaniline films for different doping degrees and different temperatures. 
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Polyaniline is a common name for a family of polymers consisting of a sequence of 
oxydized [-(C6H4)-N=(C6H4)=N-] and reduced [-(C6H4)-(NH)-(G6H4)-(NH)- 
] units. Three distinct forms have been identified and the emeraldine form, 
which contains an equal number of oxydized and reduced segments, has been 
the most frequently studied. Polyaniline was synthesized and processed as de- 
scribed elsewere [fL2[ . Flexible and free-standing films were obtained. Upon 
protonation, plunging the emeraldine sample into aqueous HC1 solutions of dif- 
ferent molarities, its resistivity changes from the characteristic of poorly con- 
ducting semiconductors (p ~ 10 10 fi.cm ) to metals (p ~ 10 _1 — 10~ 3 57. cm ). 
Circular, 12 pm in thick samples, had central gold electrodes of a diameter equal 
to 0.5 cm evaporated on both sides. Impedance measurements were carried 
out using a controlled Frequency Response Analyser coupled to a Potentiostat 
which operates in the 10 _1 to 10 6 Hz frequency range. 

Figure 3 shows measurements of real and imaginary components of ac re- 
sistivity of PANI films under conditions: a) doped in HC1 solution of 0.1 M 
at 300^ (measurement denoted by Ml); b) doped in HC1 of 0.01M at 300K 
(M 2 ); and c) doped in HC1 of 0.1M at 200-K" (M 3 ). It should be that the real 
and imaginary parts of the longitudinal resistivity p/ / are proportional to these 
measured resistivities |i3| . The sample used in the curves of Fig. 3 (b) is less 
doped than that of Fig. 3 (a), while the sample in Fig. 3 (c) is equally doped 
but it is obtained at a lower temperature. 
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4 Results and Discussion 

The curves of Fig. 3 are faithfully reproduced by our model by fixing R c and pd 
values and having a frequency-dependent interchain resistance Ri as shown in 
the insets of Fig. 4a. Fig. 4b (M2 measurement) is obtained with R c = 2.0lu 
and pd = 0.18. Finally for M3, R c = 0.5w and pd = 0.21 - the same values as 
in Mi the differences being on account of the hopping resistance Ri. At low 
frequencies Ri of M3 is ten times bigger than that of Mi (see insets of Fig. 4a 
and 4c). Then from these results we conclude that: i) the increase in the doping 
can either be simulated by diminishing R c or increasing pd and it) decreasing 
the temperature the sample resistance increases, as expected. The resistance 
R c represents the conduction mechanism along a single chain. It should depend 
only slightly on the temperature (see for example the SSH model Q), but it 
is sensitive to the doping degree. The doping parameter pd changes with the 
doping degree (as can be seen in curves 4(a) and 4(b)), but it is frequency in- 
dependent. In a highly disordered bulk structure, on the other hand, interchain 
charge transfer processes of electronic carriers can be represented by a distribu- 
tion of energy barriers W of interchain hopping (or phonon-assisted tunneling) . 
Since Ri represents, in our model, the interchain processes, its dependence on 
the frequency is explained in terms of this energy barrier distribution. A given 
waiting time r is associated with an energy barrier W by r ~ exp(W/kT). In 
the low frequency region, where / < r _1 , electron motion is hampered mainly 
by the high energy barriers (high W), which represent interchain obstacles. For 
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/ > t _1 , on the other hand, the electronic carriers become localized in small 
regions of low energy barriers, being therefore, more mobiles. This description 
agrees qualitatively with the frequency-variation of Ri shown in the insets of 
Fig. 4. To have a more formal insight on the subject, let us consider the random 
free-energy barrier model. 

This model is an elegant way to approach ac conductivity in polymeric dis- 
ordered systems. In this formalism all information about disorder is contained 
in the hopping time distribution function ip(t) . According to J. C. Dyre |L4| 
ip(t) is determined by hopping over random distributed energy. He used the 
continuous time random walk and the effective medium approximation to solve 
the model and to derive a simple expression for the complex conductivity a as 
a function of the frequency / 

4 f ]„( tmam ) 

*(f)=cHf+ jj^ i : ) } 

where C is a constant depending among other parameters on the density of 
carriers, "fmax and "f m in are two parameters related to the jump frequency. 

In our model, the Ri dependence with / was determined by adjusting the 
real and the imaginary parts of the simulated resistivity, for each frequency, to 
the experimental values. This give us an infinity of fitting parameters. But 
now we can compare our results with those predicted by Dyre's solution. The 
advantage is that Dyre's hopping resistance has only three fitting parameters for 
all frequencies. Figure 5 shows that the adjusted Ri is in reasonable agreement 
with Dyre's theory. 
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Extrapolating the dependence of Ri with / for high frequencies (unattainable 
in our equipment), as shown in the insets of Fig. 4, our simulations anticipate 
the emergence of a second plateau for the real component of the resistivity. This 
prediction is compatible with other results found for a polymeric blend [fl5|| . 

5 Summary and Conclusions 

In conclusion, we introduce a model which carefully takes into account the 
geometric aspects of the real polymeric structure. It considers the polymer 
chains being constructed from a Gaussian size distribution and in the presence of 
the excluded volume interaction. It incorporates, in a simple way, the existence 
of both intrachain and interchain charge transport mechanisms. Simulations 
based on a tranfer matrix technique and performed at the critical percolation 
threshold, reproduces the measured resistivity in PANI films in a large range of 
frequencies for different doping and temperature conditions. Estimatives of the 
critical conductivity exponent using a finite size scaling approach will be the 
object of a forthcoming paper. 
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Figure Captions 

Figure 1. The nth cell of a typical strip. Here L = 15 and we use a gaussian 
size distribution centered in Z = 6 and width Al = 6. Continuous lines represent 
the polymer chains (R c resistors) and broken lines the interchain bridges (RC 
circuits). The two thick lines are the electrodes. 

Figure 2. The longitudinal resistivity p// as a function of the frequency / 
for several values of: (a) doping parameter pd ; (b) intrachain resistance R c ; (c) 
interchain resistance Ri and fixed capacitance C = 1 nF; (d) concentration k. 

Figure 3. Dependence of the real and imaginary parts (p' and p" which are 
represented by thick and thin lines, respectively) of the measured polyanilinc 
resistivity with the frequency / for different temperatures and doping degrees. 

Figure 4. Plots of the real and imaginary parts (p'^ and p'^) of the simulated 
longitudinal resistivity (thick and thin lines, respectively) as a function of the 
frequency /. To reproduce the experimental curves, the interchain resistance 
Ri should depend on the frequency as shown in the insets. The broken lines 
corresponds to extrapolated values as explained in the text. 
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Figure 5. Dyre's hopping resistance (-Re(i)) and the adjusted Ri (inset of 
Fig. 4a) versus the frequency /. The values of the parameters C, j m in and 
j max are 1.0 10 10 , 200 and 1.0 10 9 . 
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